Fermi coordinates, simultaneity, and expanding 
space in Robertson- Walker cosmologies 



Explicit Fermi coordinates are given for geodesic observers comoving with 
the Hubble flow in expanding Robertson- Walker spacetimes, along with exact 
expressions for the metric tensors in Fermi coordinates. For the case of non 
inflationary cosmologies, it is shown that Fermi coordinate charts are global, 
and space-time is foliated by space slices of constant Fermi (proper) time 
that have finite extent. A universal upper bound for the proper radius of any 
leaf of the foliation, i.e., for the proper radius of the spatial universe at any 
fixed time of the geodesic observer, is given. A general expression is derived 
for the geometrically defined Fermi relative velocity of a test particle (e.g. 
a galaxy cluster) comoving with the Hubble flow away from the observer. 
Least upper bounds of superluminal recessional Fermi velocities are given for 
spacetimes whose scale factors follow power laws, including matter-dominated 
and radiation-dominated cosmologies. Exact expressions for the proper radius 
of any leaf of the foliation for this same class of spacetimes are given, ft is 
shown that the radii increase linearly with proper time of the observer moving 
with the Hubble flow. These results are applied to particular cosmological 
models. 
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1. Introduction 



The usual foliation of a Robertson- Walker spacetime (A4,g) by maximally sym- 
metric space slices parameterized by synchronous proper time t, deter- 
mines a notion of simultaneity and leads in a natural way to Hubble's law, 

d(t) ee v H = Hd. (1) 

Here H is the Hubble parameter, d is the proper distance on £ t from the observer 
to the test particle (such as a galaxy cluster), and the overdot on d signifies dif- 
ferentiation with respect to t. If H > and the distance d is sufficiently large, 
the Hubble speed Vh exceeds the speed of light. It is argued largely on this basis 
that in a physical sense, the universe is expanding, c.f., for example, [TJ [2j [3j H] . 

However, for a given geodesic observer with world line "/(t) following the Hubble 
flow, another geometrically natural foliation by space slices {A4 r }, previously 
considered in [5J [7J |S], is also available and leads to a different definition of 
simultaneity, and to a different notion of relative velocity. To define A4 T , let 
(p T : M -> R by, 

<Pr (p) = 9(exp~l T) P, l'{r)). (2) 

Then, 

Mr = <p-\Q). (3) 

In Eq.Q the exponential map, exp p (u), denotes the evaluation at affine param- 
eter 1 of the geodesic starting at the point p € A4, with initial derivative v. 

We refer to M T as the Fermi space slice of r -simultaneous events and to the 
observer following the timelike geodesic path j(t) as the Fermi observer. Ai T 
consists of all the spacelike geodesies orthogonal to the path of the Fermi ob- 
server at fixed proper time t = r. The restriction g T of the space-time metric g 
to A4 T makes (Ai T ,g T ) a Riemannian manifold. 

Fermi coordinates are associated to the foliation {.Mr} in a natural way. Each 
spacetime point on M. r is assigned time coordinate r, and the spatial coordi- 
nates are defined relative to a parallel transported orthonormal reference frame. 
Specifically, a Fermi coordinate system [3 [9J [TUl E] along 7 is determined by an 
orthonormal tetrad of vectors, eo(r), ei(r), e,<i{f), &z{f) parallel along 7, where 
e (r) is the four- velocity of the Fermi observer, i.e., the unit tangent vector of 
7(7-). Fermi coordinates 3 relative to this tetrad arc defined by, 

x° (exp (T) (A 3 ei (r))) = t 
x (exp 7(T) (A 3 e,(r))J - \ k , 
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where Latin indices run over 1,2, 3 (and Greek indices run over 0, 1, 2, 3). Here 
it is assumed that the X J are sufficiently small so that the exponential maps in 
Eq. (0} are defined. 

From the theory of differential equations, a solution to the geodesic equations 
depends smoothly on its initial data so it follows that Fermi coordinates are 
smooth, and it may be shown in general that there exists a neighborhood U of 
7 on which the map (a: , x , x 2 , x 3 ) : U — > K 4 is well-defined and is a difleomor- 
phism onto its image, [12] (p. 200). We refer to (x a , U) as a Fermi coordinate 
chart for 7. 

A particularly useful feature of Fermi coordinates is that the metric tensor 
expressed in these coordinates is Minkowskian to first order near the path of 
the Fermi observer, with second order corrections involving only the curvature 
tensor [TO]. General formulas in the form of Taylor expansions for coordinate 
transformations to and from more general Fermi- Walker coordinates were given 
in [T3] and exact transformation formulas for a class of spacetimes were given 
in [14l [15]. Applications of these coordinate systems are voluminous. They 
include the study of tidal dynamics, gravitational waves, statistical mechan- 
ics, and the influence of curved space-time on quantum mechanical phenomena 

m nana nana 12312111231 is]. 

In this paper, we find explicit expressions for the Fermi coordinates (a; , x , x 2 , x 3 ) 
for Robertson- Walker cosmologies when the scale factor a(t) (see Eq.(J7])) is an 
increasing function of t. We show that the Fermi chart (x a ,U) for the Fermi 
observer j(t) in non inflationary!! Robertson- Walker space-times is global, i.e., 
we may take U = M. . 

We prove that the proper radius pm t of M. T is an increasing function of r and 
is bounded above by the speed of light divided by the Hubble pararmeter at 
time r, i.e., 1/H. We show, in addition, that synchronous time t decreases to 
zero along any spacelike geodesic, Y T (p), orthogonal to the path of the Fermi 
observer at fixed proper time t — r, as the proper distance p — > Pm t - Thus, 
if the cosmological model includes a big bang, i.e., satisfies the condition that 
a(0 + ) = 0, then at any proper time r of the Fermi observer, the geodesic path 
Y T (p) in M T terminates at the big bang. In this sense, also noted in 018], the 
big bang is simulaneous with all spacetime events for this class of cosmological 
models. 

As a special case, for models with scale factors of the form a(t) = t a for some 
< a < 1, including matter-dominated and radiation-dominated cosmologies, 

3 A Robertson- Walker space-time is non inflationary if a(t) < for all t. 
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we find an exact expression for the proper radius of the Fermi space slice of 
r-simultaneous events, 



Pm 



2«r(i + i) 



(5) 



General relativity restricts the speed of a test particle to be less than the speed 
of light relative to an observer at the exact spacetime point of the test particle, 
but for test particles and observers located at different space-time points, the 
theory provides no a priori definition of relative velocity, and hence no upper 
bounds on speeds. Distant particles may have superluminal or only sub light 
speeds, depending on the coordinate system used for the calculations, and on 
the particular definition of relative velocity. 

Mitigating such ambiguities, four geometrically defined, non equivalent notions 
of relative velocity were introduced and developed in a series of papers by V. 
Bolos (c.f. [25l [26] and the references therein). Of these, the Fermi relative 
velocity is most closely analogous to the Hubble relative velocity of Eq.([T]), 
and lends itself to comparison. For the case of a test particle receding radially 
from the Fermi observer, the Fermi relative velocity may be described as follows. 

The world line of a receding test particle intersects each space slice M. T at a 
point Y T (p(r)) on the spacelike geodesic Y T (p) in M. T . The Fermi speed for such 
a particle is, 



In Eq.©, p(r) is the proper distance at proper time r from the Fermi observer 
to the test particle's position in M T . General definitions and properties of 
Fermi relative velocity for observers and test particles following arbitrary time- 
like paths are included in |26] , and Fermi relative velocities were applied in [27j 
to study the influence of a cosmological constant on receding test particles in 
the vacuum surrounding a galaxy supercluster@ 

In the present paper, we find a general expression for the Fermi speeds of test 
particles comoving with the Hubble flow relative to the Fermi observer. We 
prove that the Fermi relative velocities can exceed the speed of light in cosmo- 
logical models with a(t) — t a for some < a < 1 and find sharp upper bounds 
proportional to pm t of those speeds. We show in particular that superluminal 

4 We remark that Hubble velocity was incorrectly referred to as Fermi relative velocity in 



V F = 




(6) 
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relative velocities occur in matter dominated and radiation dominated cosmolo- 
gies, but not in the de Sitter universe, in contrast to superluminal Hubble speeds. 

The paper is organized as follows. In Sect. 2 we give basic definitions and 
solve the geodesic equations for spacelike geodesies orthogonal to the Fermi ob- 
server's path, 7(i). In Sect. 3, we construct Fermi coordinates, express the 
metric tensor in these coordinates, and prove that the Fermi chart is global for 
non inflationary cosmologies. Sect. 4 begins with the universal upper bound 
for pm t , and then gives exact expressions for Fermi relative speeds of comoving 
test particles in expanding Robertson- Walker spacetimes. Included are super- 
luminal least upper bounds for the relative Fermi speeds and exact values for 
Pm t when the scale factors follow power laws. Sect 5 illustrates and applies the 
preceding results to particular cosmological models, and Sect. 6 is devoted to 
concluding remarks and discussion of the expansion of space. 

2. Definitions and spacelike geodesies 

The Robertson- Walker metric on space-time A4 = Mk is given by the line 
element, 

ds 2 = -dt 2 + a(t) 2 (d X 2 + S 2 ( X )dn 2 ) , (7) 
where d£l 2 — d9 2 + sin 2 9 dip 2 , a(t) is the scale factor, and, 

fsin(x) if k = 1 
X if k = (8) 

sinh(x) if k = — 1. 

The coordinate t > is synchronous proper time and x, f are dimensionless. 
The values 1,0,-1 of the parameter k distinguish the three possible maximally 
symmetric space slices for constant values of t with positive, zero, and negative 
curvatures respectively. 

There is a coordinate singularity in Eq.© at x — 0, but this will not affect 
the calculations that follow. Consider the submanifold Me a ,ip = ■M-8 ,'po,k 
determined by 9 — 9q and tp = ipo- The restriction of the metric to Me , Vo is 
given by, 

ds 2 = -dt 2 + a(t) 2 d X 2 - (9) 

On Aig ^ , the coordinate x can be extended to take all real values if k = or 
— 1, and for the case that k = 1, the range of x is an interval centered at zero, 
so there is no coordinate singularity at \ — on the submanifold (see e.g., [24]). 
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Consider the observer with timelike geodesic path, ^(t) = (t, 0) in M.e , Vo . Our 
immediate aim is to find expressions for all spacelike geodesies orthogonal to 
j(t). From the Lagrangian for Eq.©, it follows that a(t) 2 dx/dp is a constant 
C along geodesies parametrized by arc length p. Since the tangent vector to the 
geodesic has unit length we get, 

dt\ 2 C 2 



The requirement that X = (t,x) is orthogonal to (1,0), the tangent vector to 
"f{t) at t = t, forces C — a(r) = oq. We will assume throughout that a(t) is an 
increasing function of t, so that by Eq. (p~0|) . t < 0, and therefore, 




a 2 (t) d X (U) 

Remark 1. The arc length parameter p for the geodesic Y(p) = (t(p),x(p)) 
with tangent vector X may be chosen so that Y(0) — (£(0),x(0)) = t(t). With 
this convention, which we assume throughout, it follows from symmetry that 
x{p) is an odd function of p. 

In light of Remark [1] there is no loss of generality in restricting our attention to 
those spacetime points with space coordinate x > corresponding to p > for 
the purpose of finding spacelike geodesies orthogonal to "/(t) with initial point 
on 7(i). 

Remark 2. A general expression for X in Cartesian coordinates was given 
in \2&l and may be deduced from Ea. \ll\ using the transformation of space 
coordinates, Sk{x) = x /(^ + \kx 2 ). 

The vector field X can be integrated to give explicit formulas for the geodesic, 
Y(p), for the special cases that a(t) = exp(H t), where H is Hubble constant, 
i.e., for the de Sitter universe [2], and for the Milne universe (see, e.g., [3]), i.e., 
for a(t) = t. To obtain integral expressions for the general case, we introduce 
a change of parameter from p to a. For the sake of clarity of exposition, we 
assume henceforth that k = or — 1 so that the range of x is unrestricted. The 
minor modifications needed to deal with the case k = 1 would include altering 
the range of a indicated in Eq. (fT2|) . but the methods are the same. Let, 



2 

a 



a(t) 



a X and a G [1, ^(t)) , (12) 



where 



/ n _ j (a(T)/a in f) 2 if ainf = Hm i _>.o+ a(t) > 
o-oo(t) = < (13) 
I oo it lim t _ ! .Q+ ait) = 0. 
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Using Eq. (fTTj) . it follows that, 



i = -VF=1, (14) 

and then differentiating Eq.lfTSf gives, 

dcr .. „a(t) , - , H 

— = a oX = 15 
dp a(i) 

where d(t) denotes da/ dt. From the chain rule, 

dp ao da dp ' 
and combining this with Eg. (1151) gives, 

dx _ aft) / 17 x 

rfcr 2a a(t)^/a - l' 

We end this section with a theorem and corollary that give explicit integral 
formulas for spacelike geodesies orthogonal to the timelike path j(t). 

Theorem 1. Leta{t) be a smooth, increasing function oft with inverse function 
b(t). Then the spacelike geodesic orthogonal to j(t) at t = t and parametrized 
by the (non-affine) parameter a is given by Y T (a) = (t(r, u),x(r, cr)) where, 



t(r,<r) = 6 (^2) (18) 



and where the overdot on b denotes differentiation. Moreover, for fixed t, the 
arc length p along Y T (a) is given by, 

, , a(r) f a • /a(r)\ 1 



Proof. Eq. (fT5)) follows immediately from Eq. (TT2"1) . To prove the other two equa- 
tions, observe that by the inverse function theorem, 

*(«(«)) = ^. (21) 

From Eq. (fl"2"]) , ait) = ao/y/a and combining this with Eq. (|2"Tj) gives, 

a(t) _ a(r) • / a(r)\ 

( j 
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Substituting Eg. ([22]) into Eqs. (fTTj) and (|T5l) and integrating yields Eqs. ([TO|) and 
1201). □ 



From Eg. ((20)) it follows that for a fixed value of r, p is a smooth, increasing 
function of a with a smooth inverse which we denote by, 

(T T {p) = <KP). ( 23 ) 

Combining Eq. (|23p with Theorem [T] immediately gives the following corollary. 

Corollary 1. Leta(t) be a smooth, increasing function oft with inverse function 
b(t). Then the spacelike geodesic orthogonal to j(t) at t = t, and parametrized 
by arc length p, is given by 

Y T (p) = (t(T,o-(p)), X (T,o-(p))). (24) 

3. Fermi coordinates 

In this section we find explicit Fermi coordinates for the timelike geodesic co- 
moving observer in Robertson- Walker cosmologies, and show that with suitable 
assumptions on the scale factor a(t), Fermi coordinates cover the entire space- 
time. We begin with some technical results needed for that purpose. 

In this section, we use the following notation, 

l7 = {(r,(T):T>0andcrG(l,<7oo(r))} (25) 

and 

C/ 1 = {(r,(7):r>0andc7e [1, ^(r))}, (26) 
where Coo (r) is given by Eq. (fT3| . Observe that U is an open subset of R 2 . 

Lemma 1. In addition to the hypotheses of Theorem^ assume that a(t) is 
unbounded and b(t) > for all t > 0. Then the map F : U\ — > (0, oo) x [0, oo) 
given by, 

F(t, a) = (t( T , a), X (t, <t)) = Y T (a), (27) 

is a bijection, and F : U — > (0, oo) x (0, oo) is a diffeomorphism. Here, the 
functions t and \ are defined by Eqs. (|18|) and (|19j) respectively. 



Proof. Let (t\,Xi) € (0, oo) x [0, oo) be arbitrary but fixed. We show that 
F(ti,o~i) — (ti,Xi) for a uniquely determined pair (ti,(7i) e U\. From Eq. (|T5)) 
we must have, 
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It remains to find t\. To that end, define the function <t(t) by, 

-ir)^(^). (29) 



a(ti 

From the hypotheses of the lemma, <t{t) is an unbounded increasing function 
of t. Then by Eq.®, 

X (r) ee X (r, a(r)) = ~ b (^ff) —^==da. (30) 



2 h 

By hypothesis, both a and b are non decreasing functions so it follows that 
x(t) is strictly increasing. A short calculation shows that x(r) has range [0, oo) 
for r > t\. Thus, by continuity there must exist a unique t\ > t\ such that 
x(t~i) — Xi- It now follows from Eq. (|2"5|) that F(tx,(Ti) — (£i,Xi) an d is a 
bijection. Now consider the restricted map, 

F : U -> (0,oo) x (0,oo). (31) 

By direct calculation, the Jacobian determinant J(r, a) is given by, 

J(Tjff) = ^f^fiM + £W r!M^,. ( 32) 



2a \ y/a J \ y/a - 1 2^ aVcr - 1 

Since a,a,b are positive and 6 > 0, it follows immediately that J(<r, r) > on 
its domain, and by the inverse function theorem F is a diffeomorphism. □ 

Let G(t, a) = (r,p(a)). Then G is a diffeomorphism with inverse, G _1 (r, p) = 
(r, ct(/3)) and non vanishing Jacobian. Using the notation of Lemma [1] define, 

H(t, x ) = GoF- 1 (t lX ). (33) 

Then H is a diffeomorphism from (0, oo) x (0, oo) onto an open subset of (0, oo) x 
(0, oo) and may be extended to a bijection with domain (0,oo) x [0, oo). We 
state this result as a corollary: 

Corollary 2. Let a(t) be a smooth, increasing, unbounded function on (0, oo) 
with inverse function b(t) satisfying b(t) > for all t > 0. Then the function 
(t, p) = H(t,x) given by j^q. (|33|) is a diffeomorphism from (0, oo) x (0, oo) onto 
an open subset of (0, oo) x (0, oo) and H may be extended to a bijection with 
domain (0,oo) x [0,oo). 

Remark 3. Using the notation of Corollary^ we may write Y T (p) = H~ 1 {t, p). 

Remark 4. The condition d(t) < for all t for a Roberts on- Walker space-time 
to be non inflationary is equivalent to b(t) > for all t. 



9 



The previous corollary guarantees that {r, p, 6, ip} are the coordinates of a 
smooth chart for the Robertson- Walker metric given by Eq.([7]). 

Theorem 2. Let a(t) be a smooth, increasing, unbounded function on (0, oo) 
with inverse function b(t) satisfying b(t) > for all t > 0. In {r, p,9,ip} coor- 
dinates the metric of Eq.^ty is given by, 

ds 2 = g TT dr 2 + dp 2 + ^Msl( X (r, a(p)))dn 2 , (34) 




2 



/ ^Lda] , (35) 



cmrf where a(p) and x(t,o~(p)) are given by Eqs. (|23|) and (|19j) . 

Proof. By the chain rule, the derivative of is given by D H -i = DfDq 1 , 



M. §1\ IdFx dF\\ (d(h dG£ 

Or dp] _ I dr da dr da 

§x dxj - [m mjldGz dG2 

><9t dp/ \ dr da I \ dr da 



(36) 



The second column of D H -i is given directly by Eqs. (|T4l) and (|T6|) . The entries 
in the first column are, 

dt_d Fl d Fl 

d^-^7 +DG - 1 ^a~ (37) 

and 

dx^dF 2 dF 2 

d^-^7 +DG - 1 ^a~' (38) 
where D 2 }_ 1 is the (2, 1) entry of Dq . Inverting the matrix Dq yields, 



a 2 Vv~T 



Using Eqs. ([27|) . (|T8j) . and (fl9|). the entries of Dp may be calculated directly and 
are given by, 

d J± = h (j^L\j^L ( 4o) 

- if °( r ) \ a ^ an 

do b \^(p))2a^ { py ^> 

dF 2 _ i r^M-n 4w =d5 (42) 



dr 2 J i \ / a\Ja — 1 
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and 



dF 2 
da 



°(t) 



(43) 

2 x /ajpj x /ajp)^l 

Using the definition of a given in Eq. (fTS)) . the first two metric components of 
Eq.Q may be expressed as functions of t and p as follows, 

a 2 (r) 



9tt 



-1 and g xx = a(t(r, p)) 



(44) 



In what follows, let x 1 , x 2 , x 3 } denote {r, p, 0, <p}. Then using Eas.lfT3]). p^|) 
and (|44l) the coefficient g pp of dp 2 in the metric tensor is given by, 



g PP — g a p - 



dx a dx@ 



1. 



(45) 



dp dp 

which may also be deduced by noting that the tangent vector 8/ dp of the 
geodesic Y T (p) has unit length (see Corollary [TJ) . Similarly, a calculation using 
Eqs.lGEB,©,©,© and 03) gives, 



dx a dx? _ d 2 (r) 



dr dr 



4a 



a 2 (r) 



a{r)\ 



2b 



i ay/ a — 1 
a[r) 



da — a 



\fa\J a — 1 



i a 2 V5^~T " o(r) Ji cr^Va^T 

2 



(«)^+K*)«_m^ 



Applying integration by parts to the integral, 



-.da, 



(46) 



(47) 



h d 3 / 2 yfa~l 

results in simplification of Eq. p5|) and yields Eg. ([331) . For the off-diagonal 
components, a calculation using Eqs. (f3"5)l . (j31?1) . (f4*2j). (H31 and (03]) results in, 



g - 1 j ) ^ a(r) 



o(r) 



(3?) 



f/<7 



cr y z Ji aVo" — 1 



a 2 \f~a~ 



1 



(48) 
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Analogous simplification using integration by parts shows that the right hand 
side of Eq. (|48|) is identically zero. □ 



The following corollary is used in the proof of Theorem [3] 

Corollary 3. With the same assumptions as in Theorem\^and using {r, p, 9, ip} 
coordinates, for fixed 8o,(fo, and r > 0, the path Y T (p) = (r, p,6o,(po) is a 
spacelike geodesic with parameter p > 0. 

Proof. We first express the geodesic Y T (p) on Me 0yipB given in Corollary [1] in 
terms of r, p coordinates. Since 

d dt d dx d 
dp dp dt dp d\ ' 

the tangent vector field along Y T (p) coincides with the tangent vector field along 
the path (r, p) with parameter p > and r fixed. Since these paths have the 
same initial data at 7(7"), they represent the same geodesic in the two respective 
coordinate systems. 

It then follows that the function Y T (p) = (t, p, 60, fo) on M is a spacelike 
geodesic. That the angular coordinates are constant follows from symmetry or 
by solving the geodesic equations directly. □ 

We refer to the coordinates, {r, p, 9, ip}, of Theorem [2] as Fermi polar coordi- 
nates. The terminology is justified by the following theorem. 

Theorem 3. Let a(t) be a smooth, increasing, unbounded function on (0,oo) 
with inverse function b(t) satisfying b(t) > for all t > 0. Define x° = t, 
x = p sin 9 cos tp, y = psin9smip, z — pcos9. Then the coordinates {r, x, y, z} 
may be extended to a chart that includes the path j(t) = (r, 0,0,0) and {d/dr, 
d/dx,d/dy,d/dz} is a parallel tetrad along 7(7"). With respect to this tetrad, 
r,x,y,z are global Fermi coordinates for the observer 7(7"). Expressed in these 
Fermi coordinates, the metric of i?g. (|34[) is given by, 



ds 2 — g TT dr 2 + dx 2 + dy 2 + dz 2 

+A fe (r, p) [(y 2 + z 2 )dx 2 + (x 2 + z 2 )dy 2 + (x 2 + y 2 )dz 2 (50) 
—xy{dxdy + dydx) — xz(dxdz + dzdx) — yz(dydz + dzdy)] , 



where g TT is given by Eq. (|35|) , p = \/ x 2 + y 2 + z 2 , and, 

p A \ k {r,p) = ^ls 2 (x(r,a(p))) p 2 . (51) 

The smooth function Afe(r, p) is a function of r and p 2 , and the notation in 
i?g. (|5ip is the same as in Theorem^ 
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Proof. We begin by showing that the indicated transformation of coordinates 
applied to Eq.(|34"|) results in Ea.([50j). Using Eq.((T6j), Eq.{(34|) may be rewritten 

as, 



ds 2 = g TT dr 2 + (dp 2 + p 2 dtt 2 ) 



a{r) S 2 k (x(T,o-( P ))) 
X 



dtt\ (52) 



where for convenience we take \ = dx/dp- Applying the change of variables in 
the statement of the theorem results in, 

ds 2 — g TT dT 2 + dx 2 + dy 2 + dz 2 

+ OfeC^P) r 2 + z , 2 + 2 + ^ 2 + 2 + y2 2 

P 

—xy(dxdy + dydx) — xz(dxdz + dzdx) — yz(dydz + dzdy)\ , 

where, 

Q h ^ P ) = a{r f^ T f p))) - P \ (54) 

From Remark [1] it follows that \ i s an even function of p, and thus Qk(j~,p) 
has a smooth extension to an even function of p. Repeated use of Eqs. (fl4]) and 
Yields, x(t,<t(0)) = 0, x(r,a(0)) = a(r)'\ x(r,a(fl)) = °: X(r,a(0)) = 
2a 2 (r) / a 3 (r) , and X (r, <r(0)) = 0, where as above the overdots signify differ- 
entiation with respect to p. Using these results, it follows by direct calculation 
that Qk(r,0) = and each of the first three derivatives of Qk with respect to 
p vanish when evaluated at (r, 0), for fc = 1,0,-1. Thus, writing Qk(r,p) as a 
Taylor polynomial in powers of p 2 , we have, 

Qk(r,p) = p i X k {T,p), (55) 

where Afc(r, p) is smooth and a function of p 2 , establishing Eq. ([50|) . which ex- 
tends by continuity to the path 7(7") = (t, 0, 0, 0), where the metric is Minkowskian. 
It now follows by calculation that all first derivatives with respect to r, x, y, or 
z of the metric tensor components vanish on 7(7"), forcing the connection coef- 
ficients also to vanish on 7(7"). Thus, each of the vectors in the tetrad {d/dr, 
d/dx,d/dy,d/dz} is parallel along 7(7"). 



Expressing the geodesic Y T (p) in Corollary [3] in terms of {r, x,y, z} gives, 

y Ap) = (t, aV, a 2 p, a 3 p), (56) 

where a 1 = sin 8 cos <po , a 2 — sin 9 sin ip Q , a 3 — cos#o- The geodesic of Eq. ([56"f 
may be extended to j(t) and is orthogonal to 7(i). It now follows from Eq.Q 
that r, x, y, z are global Fermi coordinates for the observer 7(7"). □ 
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The following definition makes some of the notation in the introduction more 
precise and will be useful in what follows. 

Definition 1. We refer to the observer following the path 7(1-) = (r, 0,0,0) 
given in the statement of Theorem as the Fermi observer. A test particle 
with fixed spatial Robertson- Walker coordinates xo,#o,¥>o and with world line 
7o( r ) = ( t j XO; $cb fo) is said to be comoving. The Fermi observer is also defined 
to be comoving. 

Remark 5. The unit tangent vector field, d/dt, of the fundamental (comov- 
ing) observers of Robertson-Walker cosmologies — i.e., the direction of time in 
Robertson- Walker coordinates — may be expressed in Fermi coordinates via, 

dt dtdr dtdp' 1 ' 

The partial derivatives in Eq.\bl} are entries of the matrix Dn(t,x) = Dq ° 
D F -i{t,x) (see Eq.§3§ and Corollary^ and may be found explicitly as integral 
expressions from -Eos. (I40[) - (|43[) . and Theorem^ 

Remark 6. The conclusions of Theorems^ and\^ continue to hold even when 
the hypothesis that b > is violated, i.e., for inflationary cosmologies, but for 
non global charts. The forms of the metric given by Eqs. (|50l) and (|34[) are valid 
in that case on some neighborhood of the Fermi observer's path ^(t). 

4. Fermi relative velocities and the proper radius of A4 T 

In this section we find a general bound for the finite proper radius of the Fermi 
space slice of r-simultaneous events, M. T , and we obtain expressions for Fermi 
velocities of comoving test particles relative to the Fermi observer. Exact results 
are given for the case that the scale factor has the form a(t) = t a for some 
< a < 1. We begin with a definition. 

Definition 2. Define the proper radius, pm t , of the Fermi space slice of r- 
simultaneous events, Ai T , by, 

^2/ b {W^w^=i d<T > (58) 



where CodY) is given by Eg. (|13p . 

The Hubble parameter, H is defined by, 

o,(t) 

H -W) (59) 
and the Hubble radius is defined to be 1/H, i.e., the speed of light divided by 
H. 
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Theorem 4. Let a(t) be a smooth, increasing, unbounded function on (0,oo) 
with inverse function b(t) satisfying b(t) > for all t > 0. Then, 

(a) at proper time r of the Fermi observer, the proper distance p to any space- 
time point along a geodesice on the space slice, M T , satisfies the inequality, 

P<Pm^<^, (60) 
and Pm t * s a monotone increasing function of time t . 

(b) synchronous time t decreases to zero along any spacelike geodesic, Y T (p), 
orthogonal to the path of the Fermi observer at fixed proper time t, as the 
proper distance p —¥ pm t > & n d t is strictly decreasing as a function of p. 

Proof. By hypothesis, b is an increasing function, so from Eg. (1201) . the proper 
distance from the observer to a spacetime point on M T corresponding to any 
parameter value a satisfies, 



2 A \Vt J a^/ 2 Vd~^l 

o(r) . . . [°° 1 



(61) 



<t) 1 2 _ 1 



2 d(r) H 

It follows immediately from Eqs. (|T3)) and (IBTj) that pm t increases with r. This 
proves part (a). To prove (b), observe that since a : (0,oo) — > (awt, oo), then 
b : (ainf, oo) — > (0, oo) and by hypothesis b is an increasing function. Thus, 
lim a _> olllf b(a) = b(a^ n{ ) — 0. Using the notation of Theorem [TJ Eq.([2"3|). and 
Corollary [TJ we have, 

hm t(r, a(p)) = lim t(r, a) = lim b = b(a+J = 0. (62) 

It follows from Eq.fTTJJ that dt/dp < except at p = 0, so t is strictly decreasing 
as a function of p. 

□ 

The following corollary is an immediate consequence of Theorem |4jo. 



15 



Corollary 4. Under the hypotheses of Theorem^ no two distinct space-time 
points are simultaneous with respect to both synchronous time t and Fermi time 
t when t = t. 

Remark 7. The Fermi hypersurfaces of r- simultaneous events, {A4 T }, are 
isotropic with respect to the Fermi observer, but they fail to be homogeneous, 
in contrast to the t- simultaneous hypersurfaces {St} for synchronous time t 
(c.f. JE}/)- This is explained by Corollary^ and Theorem^) which shows that t 
decreases with proper distance along each M T , so the intrinsic curvature is not 
constant on M T . We discuss the significance of Theorem^ in the concluding 
section. 

It is well-known that the motion of a comoving test particle follows Hubble's 
Law, 

vh(xo) = a{r)xo = Hd. (63) 

Here the Hubble speed vh(xo) = d, the Hubble parameter H is given by 
Eq. (|59[) . and d — a(r)xo is the proper distance along the spacelike path Z T (x) = 
(■^Xi^CbVo) as x varies from to xo- Both Z T (x) and the path Y T (p) given 
by Eq. (|5"6")) — described in different coodinate systems — are orthogonal to the 
path of the Fermi observer at the spacetime point 7(t), but Y T (p) is geodesic 
whereas Z T (x) is not. 

From Eqs. (fl~9| and ([20)) , it follows the the coordinate x is a smooth, increasing 
function of p along the geodesic Y T {p). We denote the inverse of that function 
(with fixed r) by p(t,x)- The Fermi speed, Uf(xo)j of the radially receding, 
comoving test particle with world line 70 (r) = (r, xo, #0, fo)i relative to the 
observer 7(7"), is given by, 

vf(xo) = -^P( T > Xo) = P- (64) 

Eq. fM)) follows from Prop. 3 in [33] and is a special case of the more generally 
defined Fermi relative velocity for test particles and observers following arbitrary 
world lines. 

Remark 8. In analogy to a well-known expression for the Hubble speed of a test 
particle with peculiar velocity, the following identity holds for the Fermi relative 
speed of a comoving particle, 

VF ( Xo )=H(r)p + a(r)^(-^r\, (65) 

as may be verified by direct calculation. The second term on the right side 
of -£7q. (|65|) is roughly analogous to the peculiar velocity in Robertson- Walker 
coordinates. 
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The following Theorem provides a general expression for the Fermi speed vp(xo)- 

Theorem 5. Leta(t) be a smooth, increasing function oft with inverse function 
b(t). The Fermi speed, i>f(xo); of the comoving test particle with world line 
7o(r), relative to the Fermi observer, is given by, 



Vf(xo) 



q(t) 



-o(r) 



f air) 



da 



da 



o(r) 



o(t)\ 



'i a 2 \/a — 1 Co 7i cr-v/cr — 1 

where ao is the unique solution to x( T ? cr o) = Xo in i?Q. (|19p . 

Proof. With xo fixed, differentiating Eq. (j!?0)) with respect to r yields, 



d d(r) 
^(r,Xo) = — 



1 cr^V 7 ^^! 



o(t) • Z' o(r) 



1 



cfcr 



2 " \V^J a^ 2 y/on~~T dr 
Now differentiating Eq.(fT9|) with respect to r gives, 

1 dao 



dr 



= 



+ ^6 

1 o'V ' — 1 ^ 



(66) 



(67) 



(68) 



Multiplying Eq. (|B5|) by a(r)/ao and solving for the last term on the right hand 
side of Eq. (|67|) gives the desired result. □ 



Corollary 5. Let a(t) be a smooth, increasing function of t with inverse bit) 
such that b > 0, then v_f(xo) is a monotone increasing function of Xo- 



Proof. Differentiating Eq. ([66|) with respect to ao gives, 

a(r)( ^ b («M) +a{T) r l &) 



2al 



V°~o - 1 



1 a\Ja — 1 



da 



(69) 



With the assumption that b > 0, we see this expression is positive on its domain. 
Since ao is an increasing function of xo, the result follows from the chain rule. □ 
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In the next corollary we consider the class of Robertson- Walker spacetimes for 
which the scale factor has the form, 

a(t) =t a < a < 1. (70) 

It is easily checked that for these models, b(t) > 0, and therefore by Theorem^ 
Fermi coordinates are global. This class of spacetimes includes the Milne uni- 
verse (a = 1), radiation-dominated universe (a = 1/2), and matter-dominated 
universe (a = 2/3) considered in the next section. The following corollary is a 
consequence of Theorem [5j 

Corollary 6. In Robertson-Walker spacetimes with a{t) = t a for < a < 1, 

the Fermi speed of the comoving test particle with world line Jo(t), relative to 
the Fermi observer is time independent and is given by, 

Mxo) = ^-( [ ! * - da+^—^ [ ; do). (71) 

The least upper bound for {i>f(xo)} * s given by, 

Iim v F ( X0 ) = fl { ^ + Jl (72) 

The right side of Eq. (|72[) is bounded above by 1/a with equality only for a = 1 . 

Proof. By Lemma [SJ vf(xo) is a strictly increasing function. It follows from 
Lemma [T] that the limit on the left hand side of Eq. (jT2l may be found by taking 
the limit of Eq. tfTTj) as ao —> a^r) = oo. An application of L'Hopital's rule 
shows that the limit of the second term in Eq. (f7"Tj) is zero. Consider the first 
integral in Eq. flTT]) • For < a < 1, 

< a (73) 



a^ +1 \Ja — 1 <j2y/a — 1 
with equality only for a = 1. Thus, 



1 <fa<2./^Hi. (74) 



a 2q 



+ i 



\Ja - 1 V c 



Therefore, 



^r(^ + i 



-do- < 2, (75) 



r(£ + l) A a^ +1 Va~T 
from which the limit and upper bound follow. □ 

The proof of the following corollary follows by direct calculation. 
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Corollary 7. In Robertson-Walker spacetimes with a(t) — t a for < a < 1, 

the Fermi speed of the comoving test particle with world line 70 (r) satisfies the 
following relationship. 

, , p a-1 f a " 1 
«f(Xo) = - + . / i . - da, 76 

where the proper distance p of the test particle from the Fermi observer is given 
by Eq.QE®. 



Corollary 8. In Robertson-Walker spacetimes with a(t) — t a for < a < 1, 
i/ie proper radius pm t °f the Fermi space slice of r- simultaneous events, Ai T , 
is a linear function of r and is given by, 

Proof. The result follows by taking the limit as o-q —> 00 of of both sides of 
Eq. (|T6"l) . applying Eq.([7^|). and observing as in the proof of Corollary [51 that 
the second term on the right side of Eg. (1751 converges to zero. □ 



Remark 9. Since the coefficient of t in Eq. \77\ is a decreasing function of 
a, the proper radii, at fixed t, of the r- simultaneous events, M T , decrease as 
functions of a in Robertson-Walker spacetimes with a(t) — t a for < a < 1. 
Conversely, pm t ~^ °o f or an V fixed r > as a — > + . 

5. Particular Cosmologies 



In this section we apply results of the previous sections to particular cosmolo- 
gies: the Milne universe, de Sitter universe, radiation-dominated universe, and 
matter-dominated universe. We include the first of these for purposes of illus- 
tration only, as the results are well known (see e.g., [3]). In the inflationary de 
Sitter universe, the expressions for the metric components in Fermi coordinates 
for a timelike geodesic observer are also known [T3] and [TS], but we show here 
that co-moving particles necessarily recede from the observer only with Fermi 
velocities less than the speed of light, in contrast to their Hubble velocities (see 
Eq . ) . To our knowledge, the results below for the radiation-dominated and 
matter-dominated universes are new. 



The Milne Universe 



The Milne Universe is a special case of a Robertson- Walker spacetime and a 
useful prototype cosmology. It is a solution to the field equations with no 
matter, radiation, or vaccum energy. For this space-time, k = — 1 and 

a(t) = t, (78) 
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and we have for the inverse function b, 



b{t) =t b = l b = 0>0. (79) 

To find Fermi coordinates for a comoving observer, we first integrate Eq. (f2"0")) 
with the result, 



p = T J?—±, (80) 



and thus, 



1 



- — 9- (81) 
lE\ 2 



It now follows from Eqs.(jT8l) and (p]9|) that 



t = — = (82) 



and 



X = \n(^ + s/a^l) = tanh" 1 = tanh -1 . (83) 

Eqs. (|82|) and (|83| are easily inverted to give, 



t =i coshx 
p =t sinhx. 



(84) 



The proper radius, Pm t i of the Fermi space slice of r-simultaneous events, 
given by Eq. ([77|) for this example is, 

which shows that the upper bound given by Theorem 0] is sharp. Theorem [3] 
guarantees that the Fermi coordinates defined by the coordinate transformation 
determined by Eqs.([8"2"f and (JS3J) a re global, and Eq. (|3"4"|) for this case reduces 
to the polar form of the Minkowski line element, 

ds 2 = -dr 2 + dp 2 + p 2 dfl 2 . (86) 



We can thus recover from Eqs. ([85l) and ([86]) the well known result that the 
Milne universe may be identified as the forward light cone in Minkowski space- 
time, foliated by negatively curved hyperboloids orthogonal to the time axis. 
The orthogonal spacelike geodesic Y T (p) expressed in Fermi (i.e., Minkowski) 
coordinates has the form given by Eq.(|56j): 
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Y T {p) = (T,a 1 p,a 2 p,a 3 p), (87) 

where (o , a 2 , a 3 ) is any unit vector in R 3 , i.e., ^(p) is a horizontal line segment 
orthogonal to the vertical time axis in Minkowski space. Now from Eqs. flTTl) and 
the Fermi speed of a comoving test particle, with fixed spatial coordinate 
X at proper time r, corresponding to parameter a, is given by, 



vf = \ = tanhx = — , (88) 

V cr r 

which by Eq. ([85|) cannot reach or exceed the speed of light. Comoving test 
particles have constant Fermi speeds proportional to their distances from the 
observer. 



Although not new for this example, Fermi coordinates for a Fermi observer in 
the Milne universe lead to an interpretation of the Milne universe that is not 
immediately available via the original Robertson- Walker coordinates. The "big 
bang" may be identified as the origin of Minkowski coordinates, and space-time 
itself may be defined as the set of all possible space-time points in Minkowski 
space that can be occupied by a test particle whose world line includes the origin 
of coordinates, i,e., the big bang. Space does not expand, rather, idealized test 
particles from an initial "explosion" merely fly apart from the Fermi observer in 
all directions with sub light Fermi velocities. Similar interpretations were given 
in0HIB]. 



The de Sitter Universe 



The line element for the de Sitter Universe with Hubble's constant Hq > is 
given by Eq.© with k = and, 

a(t) = e Hot a(t) = H e Hot , (89) 

and thus, 

6(*) = ^M*) Kt) = Tki Kt)=-T^<0- (90) 
The de Sitter universe is a vacuum solution to the field equations with cosmo- 
II 



logical constant A = 3H'i. Eas.([T5 | .([l9 l and (gOJ) yield, 



and, 



p=-^sec- 1 (^), (91) 



1 , i— ln(cos(iJ p)) 
t = T-— lnV^ = r+ v 1. 0H " , (92) 
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e- HoT tan(g p) 
X =H^ = H ■ (93) 

It follows immediately from Eq. (|9T|) that Hop < ir/2 (note that the hypothesis to 
Theorem|3]is violated here). If synchronous time t is required to be positive, then 
it follows from Eq.flHTJ together with Eqs.([l2"]) and (H3J) that y/a < exp(if i"), 
and therefore, 

e- HoT < cos(JJop)- (94) 
Thus, along the spacelike geodesies, 

x< sin^)l (95) 

Although Fermi coordinates are not global for this example, in light of Remark 
[6]we may calculate the metric coefficients in {r, p, 9, ip} coordinates by inserting 
Eqs.flHD|), (EU), and flM]) into Eqs.flM]) and (JSHJ) . The result is, 

ds 2 = - cos 2 (H p)dr 2 + dp 2 + S[n2{H 2 QP) dn 2 , (96) 

which is the same expression obtained in [T3] . A description of the way in which 
Fermi coordinates break down at the boundary of the Fermi chart in de Sitter 
space was included in |15j . It is intriguing to observe that maximal Fermi charts 
in both the Milne and de Sitter universes each occupy a single "quadrant" of 
larger embedding space-times, |28j . 

The Fermi relative speed of a comoving test particle at time r with fixed space 
coordinate X-, corresponding to the parameter cr, is given by Theorem [5] and 
reduces to, 

v F = (97) 

(7 

Combining this with Eq. ([§5)) gives, 

u e H r 

-W - l + i ke»* xr < 98 > 

Thus, the Fermi relative speed of a comoving test particle is bounded by one- 
half the speed of light for all values of Hq, \ and r. Although Hubble and Fermi 
speeds are not directly comparable, as we discuss in the concluding section, the 
above expression is strikingly different from the standard formula for the Hubble 
speed of a comoving test particle, 

v H (x) = H a e H « T X , (99) 
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which is unbounded. 



Radiation- Dominated Universe 

The radiation-dominated universe is characterized by k = and the scale factor, 

a(t) = Vt, (100) 

and thus, 

b(t) =t 2 b = 2t b = 2 > 0. (101) 
Eqs.dUl), ([111) and (J20j) yield, 

< = -, (102) 
a 

X = l^sccT 1 y/Z, (103) 

and 

p = T ^_zl +sec -i^\ . ( 10 4) 

By Theorem [5] the proper radius of the space slice of r-simultaneous events, 
M. T , is bounded by 2t. The exact value, given by Eq. (|TT|) . is, 

By Corollary [2J < and x arc smooth functions of t and p. As in Eq.(|23p we write 
a = cr T (p). Then from Eq. (|34p , the line element for the radiation-dominated 
universe in polar Fermi coordinates is, 

ds 2 = (1 + V^~T sec" 1 y/af dr 2 + dp 2 + ^ (2r sec -1 y/^f dfl 2 . (106) 

The Fermi relative speed of a comoving test particle with fixed coordinate x 
corresponding to parameter a may be calculated from Eq.(|7ip as, 



V 'a — 1 a — 1 , _ , 
v F = - + sec" 1 ^. (107) 

(7 <J 

Using (|103[) . Dp can also be expressed in terms of x as, 

Applying Corollary [SJ we find that the asymptotic limit of the Fermi relative 
speed of a comoving test particle is ir/2 times the speed of light. 

It follows from Corollary [7] or by directly comparing Eqs. (|104j ) and (|107p that, 
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p sec 

v F = 

t a 



(109) 



so that for large a, or equivalently for large proper distance p, vp ~ p/t, in anal- 
ogy to the Milne (or Minkowski) universe. However, in the Milne universe, the 
proper distance from the Fermi observer at time r is bounded by, and asymp- 
totically equal to, t. Thus, vf = p/t < 1. By contrast, the corresponding 
bound in the radiation-dominated universe, i.e., the radius of M T , is (tt/2)t so 
that vf = p/t < t/2, with asymptotic equality. From this point of view, the 
existence of superluminal Fermi velocities in the radiation-dominated universe 
may be attributed to the greater diameters of the Fermi space slices, {Ai T }, in 
comparison to the Milne space-time. 

Matter-Dominated Universe 



The final case we consider is the matter-dominated universe. For this spacetime 
k = and the scale factor is given by 



a(t) = t a , 

where a = 2/3. The inverse of a and its derivatives are given by, 
b(t) = t 3 / 2 b(t) = It 1 ' 2 b(t) = |t-V2 > 0. 

By Theorem [U 



t = 



- ^1/3 



X = 3t 



/^r(f) 2 iq(i, I; |; i) 



r(i) 



and 



r(i) 



a - 1 



3/2 



- ' 1 4'2'4'cr 



.3/4 



(110) 
(111) 

(112) 
(113) 

(114) 



where 2-F1 (•)•;■> ") i s the Gauss hypergeometric function. It follows from Theo- 
rem [3] that the metic in Fermi polar coordinates is given by, 



T 4/3 

ds 2 = g TT dT 2 + dp 2 H x(t,o-) 2 dQ 2 , 

a 



(115) 



where cr = 0>(p), as in Eq. ([23|) . is the inverse of the function given by Eq. (|114|) , 
x(r, cr) is given by Eq. (|113p . and, 



r 4/3 



r(|) 



2 J J- V 4 ' 2 ' 4 ' 



(IK 
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From Eg. (1711) . the Fermi speed of a comoving test paricle with worldline 70 is 
given by, 

3 / f a ° 1 1 f aa 1 \ 
v F ( X0 ) = - / j - da - — / a . = da . (117) 

Corollaries |5] and |S] then give the supremum of this speed and the proper radius 
of M. T as, 

hm ,M = |4^ = ^ k ^ L31103 ' ( 118 ) 
ao^too 4 r(^) T 

This shows that the matter-dominated Universe supports supcrluminal Fermi 
velocities at proper distances away from the Fermi observer, sufficiently close to 
the radius, pm t j of the Fermi space slice of r-simultaneous events. 

6. Conclusions 

Theorems [2] and [3] of this paper give transformation formulas for Fermi coor- 
dinates for observers comoving with the Hubble flow in expanding Robertson- 
Walker spacetimes, along with exact expressions for the metric tensors in those 
coordinates. We have shown that Fermi coordinates are global for non inflation- 
ary cosmologies, i.e., when the scale factor a satisfies the condition, a(t) < 0. 
Global Fermi coordinates may be useful for the purpose of studying the influ- 
ence of global expansion on local dynamics and kinematics [2]. 

Our results also apply to cosmologies that include inflationary periods, though 
in such cases the Fermi charts are local. However, if the space-time includes an 
early inflationary period, but a(t) < for all t > to, for some to, then by recali- 
brating the scale factor to a(t) — a(t + to), a global Fermi coordinate chart and 
all of our results are immediately available for the submanifold of space-time 
events with t > to- 
la Sect. 4 we found exact expressions for the Fermi relative velocities of comov- 
ing (and necessarily receding) test particles. It was shown that superluminal 
relative Fermi velocities exist, and that those velocities increase with proper 
distance from the observer. Superluminal least upper bounds were given for 
cosmologies whose scale factors follow power laws. We note that although the 
overall qualitative behavior of the relative speeds vh and vf may be compared, 
it follows from Corollary |4] that at any given proper time of the Fermi observer, 
v h and Vf measure speeds of the same comoving test particle (with fixed co- 
ordinate xo) only when it is at different spacetime points. For a comoving test 
particle at a given spacetime point, Vh and vp give the particle's relative speeds 
at different times of the Fermi observer. 
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The existence of superluminal relative velocities bears on the question of whether 
space is expanding, c.f. [IJ [2j [3j HJ [27] and the numerous references in those pa- 
pers. On this matter, our results may be contrasted with arguments given in 
pjj. In that paper, the coordinate transformation for the Milne universe, re- 
peated in our Eqs. (j55|) and ||81j| , was used to compare the Hubble and Fermi 
(or Minkowski) relative speeds vh and vf of comoving test particles. The in- 
congruity of superluminal Hubble speeds and necessarily subluminal Minkowski 
speeds in the Milne universe was discussed. It was argued that the analogous 
qualitative difference would also occur for cosmologies that include matter or 
radiation, through a comparison of Hubble speeds to speeds defined via coordi- 
nates with a "rigid" radial coordinate^ However, if the latter class of coordi- 
nates includes Fermi coordinates — the coordinates used to deduce that conclu- 
sion for the Milne universe — our results in Sect. 5 for the radiation-dominated 
and matter-dominated cosmological models do not support that conjecture. 

For cosmological models with a scale factor of the form a(t) = t a for < a < 1, 
the existence of superluminal relative velocities of comoving particles may be 
understood in terms of the geometry of the simultaneous space slices, {M T }. 
In Sect. 5, it was shown through the use of specific examples that superluminal 
relative Fermi velocities exist provided "there is enough space" in the sense that 
the proper radius pm t of M T satisfies the condition pm t > t. 

In Sect. 5 we also showed that the relative Fermi speeds of comoving particles 
are bounded by one-half the speed of light in the de Sitter universe, a space-time 
considered to be "expanding." By way of contrast, superluminal relative Fermi 
velocities were proven to exist in |27j in the static Schwarzschild space-time 
(with interior and exterior metric joined at the boundary of the interior fluid), 
a space-time that is not usually regarded as "expanding." Thus, it may be 
argued that existence of superluminal relative velocities, in general, is not the 
appropriate criterion for the purpose of defining what is meant by the expansion 
of space. 

Does space expand? An affirmative answer may be given for the non inflationary 
Robertson- Walker cosmologies studied in this paper, in the following sense. For 
any comoving geodesic observer, the Fermi space slices of r-simultaneous events, 
{.A/fr}, that foliate the space-time have finite proper diameters that are increas- 
ing functions of the observer's proper time. This is Theorem [4^,. Theorem [4}} 
explains how it is possible for the space slices to have only finite extent. What 
stops these hypersurfaces of constant Fermi time from continuing beyond their 
proper diameters? The theorem makes precise the way in which all space-time 
events are simultaneous to synchronous time t = 0, the big bang in cosmological 

5 c.f. Sect. V of 0], p. 63. 
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models admitting an initial singularity. 
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